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ABSTRACT 

We explore the sensitivity of massive stars to neutrino magnetic moments. We find that the additional 
cooling due to the neutrino magnetic moments bring about qualitative changes to the structure and 
evolution of stars in the mass window 7M Q < M < 18 M Q , rather than simply changing the time 
scales for the burning. We describe some of the consequences of this modified evolution: the shifts in 
the threshold masses for creating core-collapse supernovae and oxygen-neon-magnesium white dwarfs 
and the appearance of a new type of supernova in which a partial carbon-oxygen core explodes within 
a massive star. The resulting sensitivity to the magnetic moment is at the level of (2 — 4) x 10~ n ^b- 
Subject headings: neutrinos — stars: supergiants — stars: evolution — stars: interiors 



1. INTRODUCTION 

Cooling by neutrino emission is known to be impor- 
tant in a variety of stellar systems, from low-mass red 
giants (RG) and horizontal branch (HB) stars to white 
dwarfs (WD), neutron stars and core-collapse supernovae 
(CCSN). Although, compared with photons, neutrinos 
interact extremely weakly, once produced they easily 
escape from stellar interiors carrying away energy that 
would otherwise take much longer to be transported to 
the surface by radiation or convection. The resulting en- 
ergy sink in the center of the star can dictate the star's 
rate of nuclear burning, structure and evolution, and ul- 
timately how it ends its life. If a new interaction modifies 
the rate at which neutrinos are produced, the evolution 
of the star could change. Stars can therefore be used as 
laborator ies to search for non-s tandard properties of the 
neutrino (|Bernstein et al.lll963f) . 

One such new interaction that will be considered in 
this paper are the neutrino (transition) magnetic mo- 
ments. The magnetic moments couple neutrinos to pho- 
ton through the effective Lagrangian term 



L= ~^fJ-l J TpiOapTpj F' 



a/3 



(1) 



where tp is the neutrino field, F the electromagnetic field 
tensor, a, (3 are Lorentz indices, i, j are the flavor in- 
dices. This interaction is non-zero, although extremely 
small, already in the Standard Model (SM) with non- 
zero neutrino mass. In some extensions of the SM it 
is allowed to be as large as the present direct labora - 
tory bound of \i v < 10~ 10 \ib (jDaraktchieva et aT1l2005l ). 
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where {1b = e/2m is the Bohr magneton (see Sect. 12.1) ). 
Such new physics could reside just above the electroweak 
scale of 100 GeV. 

The effects of the magnetic moment interaction on the 
evolution of low-mass stars have been extensivel y stud- 
ied, as will be rev iewed in Sect. [2] (see also, e.g.. iRaffeltl 
119961 11999L 120001 for a review). In fact, these stars are 
thought to give the best bound on the neutrino magnetic 
moment, about an order of magnitude more restrictive 
than the direct laboratory bound. The effect of the mag- 
netic moment on the evolution of massive (M > 7M ) 
stars, however, has not been established. This paper is 
an exploratory investigation designed to fill this gap. Our 
goal is to answer the following questions: 

• Does the presence of a non- vanishing neutrino mag- 
netic moment influence in any way the behavior of 
a massive star? 

• If so, for which range of stellar masses? 

• What are the induced effects? 

• What is the sensitivity of massive stars to the size 
of the neutrino magnetic moment? 

The paper is organized as follows: In Sect. [2j we dis- 
cuss the neutrino magnetic moments from the theoretical 
and phcnomcnological points of view. In particular, we 
review the present bounds from experiments and astro- 
physical considerations. In Sect. [3l we review the stan- 
dard evolution of massive stars and consider how might 
change with the addition of the extra cooling induced by 
a non-zero fi u . In Sect. |4j we present the results of our 
numerical modeling. Finally in Sect. we summarize 
our results and discuss possible future directions. 
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2. NEUTRINO MAGNETIC MOMENT 

2.1. Theoretical considerations 

The search for the neutrino magnetic moment 
go es back to the neu trino discovery experiment 



bv lCowan fc Reinesl (|1957f) . who put the first bound on 
the magnetic moment, \x v < 10 -9 /xb 

On the theoretical front, it was realized many decades 
ago that a neutrino with a non-zero mass should 
have a non-vanishing magnet ic moment even with the 
Standard Model inte r actions dFuiikawa fc Shrock 1 119801 : 
iLee fc Shrockl [l977t iMarciano fc Sandal I197tTT ^ Due 
to the specific nature of the SM interactions, however, 
namely, coupling of W to left-handed currents only, the 
induced value of the magnetic moment turns out to be 
very small, fj, v / /jb - 3 x 10~ 19 x (m„/leV). In exten- 
sion of the SM, this limitation need not apply and much 
larger values of \i v can be obtained. 

There are several classes of models beyond the SM pre- 
dicting a large magnetic moment (as large as /i^/^B ~ 
10~ 10 ). These scenarios typically manage to simulta- 
neously have a small neutrino mass and a large mag- 
netic moment by one of the two following mechanisms: 
(i) requiring that the new physics dynamics obeys some 



(approximate) symmetry t hat forces 



while 



allowing non-vanishing ji„ (IBabu fc Mohapatra 
Barbieri fc Mohapatra I [1989: Georgi fc Randall 





199C 




1990 



Grimus fc Neufeld Ill99ll : IVoloshin Ill988f k or by (ii) en- 
gineer ing a spin suppr ession mechanism that keeps m„ 
small (jBarr et al.lll990f ). 

It is important to stress that the value of the mag- 
netic moment in SM extensions cannot be arbitrarily 
large. It can be shown in a model-independent way that 
a large value of \x v would induce, via electroweak ra- 
diative corrections, a neutrino mass above the current 
experimental limit (conservatively, ~ 1 eV) . Quite in- 
terestingly, the resulting "naturalness" upper limits are 
substantially different for Dirac and Majorana neutri- 
nos, due to the different nature and breaking of the ap- 
proximate symmetries imposed to simultaneously ensure 
m u ~ and fx v ^ 0. For the Dirac case, a general ef- 
fective theory analysis leads to (jBarbieri fc G. Fiorentinil 
fl988tlBell et ani2005h : 



^ <3x HT 15 
Mb 



lcV 



lTcVV 



A 



(2) 



where <5m„ is the induced mass term and A is the (un- 
known) new physics mass scale at which the magnetic 
moment originates. For the Majorana transition mag- 
netic moments one finds (|Bell et al. 2006; Davidson et al.l 
120051) : 



MB 



< 4 x 10" 



[Si 



1 a d 



leV 



lTeV V 



1,1 T 



(3) 



where (5m^) a p is the induced contribution to the neu- 
trino mass matrix and m a ,mp represent the charged 
lepton masses for the flavors a, (3. Therefore, the natural- 
ness limit on the magnetic moment of Majorana neutri- 
nos is considerably weaker. Note that in both cases there 
is a strong quadratic dependence on the new physics 
scale A, for which we have taken the reference value of 
A = 1 TeV. 



In summary, evidence of [i v > 10~ 19 {1b would indicate 
the need for a radical extension of the weak sector of 
the SM. Moreover, evidence of a magnetic moment in 
the window 10~ 15 /iB ^ M„ 10~ 10 /is would provide 
strong evidence for the Majorana nature of neutrinos. 
Therefore, astrophysical probes of the neutrino magnetic 
moment with sensitivity right below /i„ ~ 10~ 10 /zb are 
of considerable theoretical value. 

2.2. Current Bounds 

The first experimental bound on the magnetic moment 
of the electron neutrino, [i v < 10~ 9 /^b, was obtained 
bv lCowan fc Reinesl (|1957f ). by analyzing electron recoil 
spectra in (anti)-neutrino electron scattering. The same 
technique has led, over time, to more stringent bounds , 
U V < (5 -10) x 10"" Mb, at 90 % C.L. dBeda et alJl2"007l : 
iDaraktchieva et aT1l2005h . A similar bound {\x v < 5.4 x 
10~ n ^B, at 90% C.L. ) was recently confirmed by the 
Borexino collaboration (Galbiati 2008]) from the analysis 
of solar neutrinos. 

Another set of constraints comes from the analy- 
sis of en ergy loss in stars. This method was pio- 
neered by iBernstein et alJ ()1963f ) who considered the ef- 
fect of the addition cooling due to the neutrino mag- 
netic moment on the lifetime of the Sun and esti- 
mated that /i„ < 10~ 10 /iB- A recent bound from 
the Sun is \i v < 4 x 10- 10 /! B (Raffclt 1999). A 
more stringent bound, [i v < 3 x 10~ 12 /is, comes 
from the analysis of the RG branch, where the ex- 
tra cooling could delay the onset of the helium flash 
(ICastellani fc deglTnnocentil 119931: ICatelan et all 119961 : 
lHaft et all 119941: iRaffeltl 119901: iRaffelt fc Weiss! I1992T . 
A comparab le bound also follows from the analys is of 
WD cooling (jBlinnikov fc Dunina-Barkovs kava 1994]). In 
contrast, the effect of a non-zero neutrino magnetic 
moment on th e cooli ng of a neutron star, studied by 
llwamoto et alJ (|1995| ). leads to the considerably less 
stringent constraint fi v < 5 x 10 _7 ^b, though a value 
as small as ~ 10~ 10 /ib would change the cooling of the 
crust and could be observable in the future. Finally, HB 
stars are also quite sensitive to a non-vanishing mag- 
netic moment. A value of [i v of about 1.5 x 10 _11 /ib 
would shorten the lifetime of the He-burning stage of 
those stars, and consequently change the number ratio of 
HB stars versus RG s tars in globular clu s ters above th e 
observational bound (|Raffelt et al.lll989l ; IRaffeltl [l996h . 
An important feature of these constraints is that they 
do not depend on the nature - Dirac or Majorana - or 
flavor of the neutrinos. 

For Dirac neutrinos with a non-vanishing mag- 
netic moment, the left handed states could be trans- 
formed into right handed states in an external elec- 
tromagnetic field. For a supernova (SN), where left 
handed neutrinos are trapped and right handed neu- 
trinos are not, this would mean a very efficient cool- 
ing mechanis m. This leads to the bound ji v < 
3 x I P" 12 «r (lAvala et al.l 11999): [Barbieri fc Mohapatral 



Il988t ILattimer fc CoopersteinI 119881 ). A comparable 
bound comes from the observation that some of the right 
handed neutrinos could transform back into the active 
ones in the galactic magnetic field, and should have been 
observed by the detec tors which me asured the neutrino 
signal from SN 1987A (|Notzoldlll988l) . 
Similarly, the analysis of big bang nucleosynthesis puts 
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Fig. 1. — Standard model evolution of the central temperature 
and density for a 25 Mq, 15 Mq and 10 Mg star. The different 
burning stages are indicated by the corresponding chemical symbol. 
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Fig. 2. — The contours of the ratio P c _|_ 7 /P e ; g of the combined 
electron and photon pressure P c + 7 to the corresponding electron 
ideal gas pressure, P e ,i.g. = n e kT, as a function of temperature 
and density of the medium. 



limits on the magnetic moment for Dirac neutrinos. A 
large magnetic moment would bring right handed neutri- 
nos in thermal equilibrium before the cosmological nu- 
cleosynthesis, spoiling the prediction of light elements 
abundances. The current b ound is approximately [i v < 
10~ n fiB (see[Dolgov 20021 and references therein). 

Lastly, a Majorana neutrino with a non-vanishing 
magnetic moment could be transformed into an anti- 
neutrino o f a different flavor in an external electromag- 
netic field. iMiranda et al.l (|2004[ ) inferred the bound of a 
few 10 _12 /zb from the lack of the observation by Kam- 
LAND of t he antineutrino fl ux from the Sun. Using the 
same data, iFriedlandl {2005), however, argued that the 
bound is actually an order of magnitude weaker. 

To conclude this Section, we point out that the bounds 
discussed above are sensitive to different effective mag- 
netic moments, namely to different combinations of the 
elements of the magnetic moment matrix fJ,"@. The ter- 
restrial experiments using reactor sources (electron anti- 
neutrinos) are sensitive to = J2 a \nt a \ 2 - On the 
other hand, the bounds based on solar neutrinos analyses 
are sensitive to a weighted average of {n v ) 2 e , (a*i/)^> (a^) 2 -) 
with roughly equal weights determined by the oscillation 
probabilities. Finally, the bounds based on energy loss 
in stars are sensitive to = ^ Q o |/i"^| 2 . 

3. MASSIVE STARS AND NON-STANDARD COOLING: 
ANALYTICAL CONSIDERATIONS 

3.1. Review of the standard evolution and cooling 

In this Subsection, for the convenience of the reader, 
we briefly review the main features of the standard evo- 
lution of massive stars relevant to our subsequent discus- 
sion of the effects of the neutrino magnetic moment. We 
collect some of the technical details in Appendix A. For a 
standard introduction t o the evolution of massive stars, 
see, e .g., lArnettl (|1996l ); iClavtonl (|1983f ): iWooslev et ail 
(|200l . 

We consider here stars in the mass range of 7M < 
M < 25 M Q . The basic information about the physical 
conditions in the core of these stars can be obtained from 
simple considerations of hydrostatic equilibrium and bal- 
ance between energy generation and loss. Assuming for 



the moment that the interiors of these stars can be de- 
scribed as ideal, nondegencrate gas, from the consider- 
ation of hydrostatic equilibrium one can find a simple 
relation among the mass M, central temperatu re T c and 
central density p c (|Arnettlll996l : IClavtonlll983h : 



T 

Ik 



4.6 x 10V 



M 



2/3 



1 g cm 



1/3 



(4) 



where fx = (£\ Yi + Y c )~\ with Y e the mean number of 
electrons per baryon, and Yi = Xi/Ai with Xi and Ai 
the mass fraction and atomic number of species i respec- 
tively. 

The curves in Fig. [T] represent the evolution of the cen- 
tral temperature and density for stars of three different 
masses, as found with our numerical code. Once one fuel 
is exhausted, the star follows a period of compression, 
which lasts until the temperature and density satisfy the 
required condition for the ignition of another nuclear re- 
action. Wc can see that these stars indeed approximately 
follow Eq. ([4]): in particular, the more massive ones have 
higher central temperatures. This simple fact will have 
important implications later. 

The approximation of ideal nondegenerate gas is valid 
for a relatively limited set of densities and temperatures, 
as seen in Fig. [2] which shows the ratio of the combined 
electron and photon pressure P e +-P ra d = Pc+j to the cor- 



responding electron ideal gas pressure, -P e ,i. ( 



-kT. At 



high temperatures and low densities the pressure is dom- 
inated by radiation (photons and relativistic electron- 
positron plasma). At lower temperatures and high den- 
sities, the pressure is set by electron degeneracy. 

Comparing Figs. [1] and [21 we see that the evolutionary 
trajectory of the 25 M Q star indeed lies almost entirely in 
the region where the ideal gas equation of state applies, 
except at the very last stages of the evolution which are 
on the edge of being degenerate. The lower mass stars 
considered in Fig. [TJ having lower central temperatures, 
are more subject to the effects of degen eracy which set 
in after the start of the carbon burning (|Arnettlll996T ). 

Two important observations need to be made at this 
point. First, so long as there is steady nuclear burning 
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at the center, the star has to be outside, or at least on 
the edge, of the degenerate region. Indeed, a strongly 
degenerate system does not have the necessary negative 
feedback mechanism to regulate the rate of burning: en- 
ergy production does not lead to expansion and cooling 
of the central region if the pressure is independent of 
temperature, a s is the case for strongly degenerate sys- 
tems (see, e.g., lRaffeltlfl996h . The detour of the 10 M Q 
star through the region of high degeneracy seen in Fig. Q] 
indicates that the nuclear burning at this stage happen s 
in a shell rather than in the core dWooslev et al.ll200l) . 
When the burning reaches the center, the temperature 
jumps and the feedback is enabled. As we will see, this 
behavior will be even more pronounced in the presence 
of the neutrino magnetic moment. 

Second, there is an important physical distinction be- 
tween the regimes of non-relativistic and relativistic de- 
generacy. The former is characterized by an adiabatic 
exponent 7 ac i = 5/3 (polytropic index n = 3/2)^ while 
the latter has 7 a d = 4/3 (n = 3). A non-relativistic de- 
generate gas responds more strongly to compression than 
a relativistic one. In fact, 7 a( j = 4/3 corresponds to the 
edge of stability: if a sufficient fraction of the stellar core 
has this adi abatic exponent, the star b ecomes unstable 
to collapse (|Zeldovich fc Novikov Ifl97ll ). In the case of 
white dwarfs, this leads to the well-known upper limit on 
the mass, the Chandrasekhar mass. In our case of extra 
cooling, this will lead to the appearance of a new type of 
supernova. 

We now turn to the mechanisms of energy loss. Energy 
generated by nuclear reactions is lost either through the 
surface of the star as photon radiation, or through the 
volume via neutrino emission. Photon energy loss per 
unit mass is about 10 4 — 1 5 ergg" 1 s -1 for stellar masses 
considered here (see, e.g.. lArnettll!996l § 6.4). Neutrino 
cooling becomes important when the emission rate per 
unit mass exceeds this value. 

At the temperatures and densities of interest for our 
discussion, neutrinos are dominantly produced by the fol- 
lowing processes (see lBeaudet et alll967t lDicuslll972l and 
reference therein): 

1. photo process, ^ e~ — > evv; 

2. pair process, e + e~ — > vv; 

3. plasma process, 7 — > vv, 
also known as plasmon decay; 

4. bremsstrahlung, e~(Ze) — > (Ze) e~ vv. 

The diagrams for the first three of these processes are 
depicted in Fig. 0] 

Which of these processes dominates depends on the 
density and temperature, as seen in Fig. [3^. The con- 
tour lines show the rate of energy loss, and are computed 
us ing the exp r ession s for the neutrino cooling rates given 
in lltoh et all ([1996) . These rates are used as a starting 
point in our numerical analyses, to compute the reference 
models with standard particle physics. 

Superimposing Fig. [1] to Fig. [3^,, we see that in the 
early burning stages (H- and He-burning), the neutrino 

t We use the standard definition P = Kp" 1 ^ = Kp 1+1 / n relating 
density p and pressure P. 



production is dominated by the photo process. The en- 
ergy loss (see the Eq. (|A1|) or contour lines in Fig. [3k.), 
however, is still small compared to the stellar luminos- 
ity, 10 4 — 10 5 ergg -1 s _1 . The neutrino energy loss starts 
dominating over the radiation only after central helium 
depletion, for temperatures above ~ 5 x 10 8 K. When 
carbon is ignited, T ~ 7 x 10 8 K, the star is already los- 
ing most of its energy through neutrinos. For the more 
massive stars (see the curve corresponding to 25 M Q and 
15 Mq in Fig. [I}, neutrinos will be mainly produced by 
the pair process. For lower mass stars, however, there 
will be a competition between the pair and plasma pro- 
cesses. As we see the bremsstrahlung process is never 
important for a massive star. 

Since the neutrino emission rate increases rapidly with 
temperature (contours lines in Fig. GJt,), the fuel is con- 
sumed faster in the late stages of stellar evolution. For 
this reason, these later stages are characterized by pro- 
gressively shorter time scales. In fact, a massive star 
spends about 90 % of its life burning H, and almost all of 
the remaining 10 % burning He. The time scale for the 
other burning stages is thousands of years for C. year s 
for Ne and O, and only days for Si (jWooslev et al.ll200"2T ). 

In our numerical model (see details in Sect. QJ, stars 
whose initial mass is above ~ 9 M Q can ignite all burn- 
ing phases (H, He, C, Ne, O, and Si burning) 1 . They end 
up with a typical "onion structure" with iron in the cen- 
ter and several layers around, in which fusion of lighter 
elements occurs. No energy can be gained by fusion reac- 
tions in the iron core, and once its mass exceeds ~ 1.4 M Q 
the star collapses producing a CCSN. 

Stars less massive than about 9M Q have a different 
destiny. If the initial mass is less than about 7.5 M Q 
they never reach the stage of carbon-burning and end 
their life as carbon-oxygen white dwarfs (CO-WDs, see 
the bottom row of Fig. [7]). If their initial mass is in 
the range 7.5 < M/Mq < 9, they do ignite carbon but 
not neon, and end their lives as oxygen-neon-magnesium 
white dwarfs (ONeMg WD). These numbers are impor- 
tant to us because they will be shifted by the presence of 
the additional cooling. Moreover, a new evolution possi- 
bility will be opened up. 

We also note that the stars that we just classified 
as making ONeMg WDs could yet have another fate. 
They are close to the Chandrasekhar mass, and, after 
a typical "deep third dredge-up" that mixes the enve- 
lope with helium layer of the core, the core may con- 
tinue to grow by hydrogen/helium shell burning. Even- 
tually, collapse of the core may occur due to electron 
capture and a su pernova can result (electron capture su- 
pernova, ECSN. lMivaii et aIJll980h . Such a supernova 
typically would produce only very little 56 Ni. The de- 
tails of whether the core can grow highly depends on 
the dredge-up (mixing) after and between the thermal 
pulses of the hydrogen/helium shell burning, competing 
with the mass loss of star that terminates core growth 
when the envelope has been lost. The h istorical devel- 
opmen ts on this topic are summarized in IWooslev et ail 
(120021) . while some recent de velopments are discussed in, 
e.g.. lPoelarends et al.l ((2008). Due to the uncertainty of 
the mass range of ECSN, we will consider them as part of 

For a discussio n of how this threshol d value varies between 
different codes, sec (Poelarcnds et al. 20QD 
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Fig. 3. — Neutrino energy loss. From above Panel a (top): 
Standard energy loss, [i u = 0; Panel b (middle): Energy 
loss for ^t„ = 10 — 11 ^b; Panel c (bottom): Energy loss for 
li v = 5 X lO -11 /te- The different colors show the regions where 
the relevant process is 90 % or more of the total neutrino cool- 
ing: Turquoise (top) neutrino pair production e + e~ — » vv\ Or- 
ange (middle left) neutrino photo production 7e~ — > evv\ Green 
(middle low) neutrino production by plasmon decay j — ► vv; Pur- 
ple (bottom right corner) bremsstrahlung production e~(Ze) — > 
(Ze) e~ uV. For comparison, sec Fig. C.4 in Raffclt (1996). 
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Fig. 4. — Feynman diagrams for the standard model neutrino 
production processes relevant for our discussion. From top to bot- 
tom: photo process, -ye~ — > euV; pair process, e + e~ — + vv; 
plasma process, 7 — » vv. 



our ONeMg core "bin" and just note that they comprise 
some fraction of the upper end of the ONeMg WD mass 
range. 

Also note that we only consider single star evolution 
here. Special evolution paths that may occur in binary 
star systems, like, e.g., Type I SNe, are not considered 
in this paper. 

3.2. Effects of the (i^-induced extra cooling 

If neutrinos have a non- vanishing magnetic moment, 
additional processes shown in Fig. [5] must be considered. 
It turns out that the main effect of the neutrino magnetic 
moment is to increase the plasmon decay rate. For ex- 
ample, for fJL u /^B = 1CP 10 , density p ~ 10 4 gcm -3 and 
temperature T = 10 s K, the non-standard plasma pro- 
cess ener gy l oss is about 10 3 larger than the standard one 
(see Eq. IB1[) . On the other hand, for the same value of 
\i v the non-standard pair production rate never exceeds 
30% of the standard one (see Fig. IB9j) . Details of the 
analysis of the modified energy loss rates are provided in 
Appendix B. 

Graphically, the effects of the magnetic moment are 
shown in Panels (b) and (c) of Fig. [3] We see that the 
plasmon dominated region enlarges considerably. Com- 
paring Fig. [3] with Fig. [T] we see that the 25 M Q star 
"misses" the plasmon dominated region, while the 15 Mq 
and 10 Mq stars, having lower central temperatures, en- 
ter it. Therefore, for them the effects of the anomalous 
cooling should be important. 

The above argument does not take into account pos- 
sible feedback effects of the extra cooling on the evo- 
lutionary trajectory. Indeed, such feedback might be 
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Fig. 5. — Feynman diagrams for the electromagnetic contribution 
to the neutrino production processes relevant for this paper. From 
top to bottom: photo process, 7C — » ev V; pair process, e+ e~ — » 
vv\ plasma process, 7 — > vv. 
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Fig. 6. — The evolution of the central temperature and density 
with the extra cooling due to neutrino magnetic moment (p v = 
0.5 X 10 — 10 /xb)- Four representative stellar masses (11.5 Mq, 
15M Q , 17M , and 20 Mq) were chosen. The 11.5 Mq star un- 
dergoes a thermonuclear CO explosion. 



role in our analysis. 2 For t he plasma p roces s we used 
the fitting function given in lHaft et all (|1994l ) which is 
shown to be better than 5 % i n all the region o f interest 
for our purpose (see Fig. 8e of lHaft et al.Hl99l . 

We used the modified code to compute a two dimen- 
sional grid of data points, varying the mass of the star, 
4 Mq < M < 25 Mq, and the value of the neutrino mag- 
netic moment, < fx v 
discussed below. 



< 0.5 x !CT lu /iB. The results are 



anticipated on general grounds: extra cooling, if suffi- 
ciently strong, can lead to the decrease of temperature 
which would in turn push the star deeper in the plas- 
mon dominated degenerate region. Once the effects of 
degeneracy become important, cooling does not change 
the pressure consistently and so there is not m uch com- 
pressional heating as a reaction (|Arnettj["l996t ). Thus, a 
qualitative change in the evolutionary trajectory might 
be expected. 

We will now see how these qualitative expectations are 
borne out by the details numerical simulations. 

4. MASSIVE STARS AND NON-STANDARD COOLING: 
RESULTS OF NUMERICAL MODELING 

4.1. Numerical modeling 

Our numerical modeling is based on KEPLER, a 
one-di mensional implicit h ydrodynamics stellar evolution 
code ( We aver et alJfl97ct ). The physical ingredients of 
this code - such as the treatment of convection, opaci- 
ties, etc - in the case of the st andard physics are well 
documented in the literature (iHeeer fe Woosievl 120081 : 
IWooslev et aT1l200l IWooslev fc Hegerll2007l) . 

In order to simulate the effects of non-standard cooling 
in the evolution of massive stars, we modified the neu- 
trino cooling routine of the code. Specifically, we com- 
puted the energy loss induced by the electromagnetic pair 
production (see appendix) and derived an interpolating 
function which is accurate at the level of 5 % in all the 
region of interest. The bremsstrahlung and photo pro- 
cesses were left unchanged since they do not play any 



4.2. Results: Overview 

Figure [6] displays the evolution of the central tempera- 
ture and density for stars of four representative masses: 
11.5M , 15M , 17M Q , and 20 M Q . These simulations 
include our largest extra cooling of [i v = 0.5 x 10 _10 /xb- 

As expected, we see that the evolution of the most 
massive of these stars, 20 Mq, shows no visible departure 
from the standard trajectory (Fig. Q}. The lower mass 
stars, in contrast, deviate from their standard trajecto- 
ries into the region of lower temperature. The curves 
start to bend between the He- and the C-burning stages, 
that is as soon as the neutrino cooling becomes more 
important than the radiation energy loss. The deviation 
becomes more pronounced the lower the mass of the star, 
again as anticipated. 

The physics behind this phenomenon is easy to un- 
derstand. As the CO core contracts following He burn- 
ing, the heat that would otherwise increase the temper- 
ature to the point of carbon ignition is instead lost due 
to the additional cooling. The core is unable to stay 
on the "ideal gas" trajectory and instead collapses to 
a degenerate configuration. The center of core is too 
cold to ignite carbon burning, which is instead ignited 
off-center. The center becomes a relatively "cold" white 
dwarf (T ~ (3 — 5) x 10 8 K) surrounded by a hotter burn- 
ing shell (T ~ 1 x 10 9 K). Notice that while off-center 
ignition is observed in the standard case for lower mass 

2 A simulation with the bremsstrahlung and photo induced cool- 
ing doubled showed, as expected, no sensitive changes in the evo- 
lution of the star. 
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stars (~ 10 M©, see Fig. [IJ, here it is seen for stars of 
15 M Q . 

What happens next depends on the mass of the star. In 
the case of the 15 M© star, the burning eventually reaches 
the center. At this point, the corresponding curve in 
Fig. [6] abruptly jumps back up to the standard "ideal 
gas" trajectory. The silicon burning stage in the cen- 
ter ensues and the star eventually undergoes iron core 
collapse. The remarkable detour into the degenerate re- 
gion results in a somewhat different composition of the 
supernova progenitor (see Sect. 14. 5j ). 

A dramatically different fate awaits stars of lower 
mass, as represented by the 11.5 M curve. In this 
case, the detour in the degenerate region is more pro- 
nounced, with the central temperature dropping down 
to T ~ (1 — 2) x 10 8 K. The shell burning of CO ex- 
tinguishes before reaching the center. As a result of the 
outer shell burning, the relativistic degenerate CO core 
builds up until it reaches the Chandrasekhar limit . The 
core then undergoes a thermonuclear runaway of carbon 
burning, consuming the CO fuel (CO explosion). This is 
indicated as indicated by the vertical line in the Fig. [6l 
In effect, the extra cooling leads to new kind of super- 
nova - a SN Type la-like thermonuclear explosion inside 
a massive star - which perhaps could be viewed as new 
"Type 1.7". Our simulation stops at this point. 

The fates of the stars of different masses as a function 
of the neutrino magnetic moment are shown in Fig. [7] We 
see that neutrino magnetic moment appreciably shifts 
the mass thresholds between the different outcomes. In 
the presence of the extra cooling, more and more mas- 
sive stars are unable to ignite carbon burning. The CO 
window expands as a result and the ONeMg window 
shrinks. The CO explosion window first opens up at 
{i v = 0.2 x 10~ 1c Vb and for \l v = 0.5 x 10 _10 /iB covers 
star with masses 9.8 M Q < M < 11.5 M©. A complete 
overview of all the models computed is given in Table [Cll 
of Appendix [Cl 

4.3. Changes of the WD/CCSN Thresholds 

As mentioned earlier, one of the effects of the extra 
cooling is to make more difficult the ignition of various re- 
actions, in particular carbon, oxygen and Neon. To burn 
those nuclei now requires a larger stellar mass, mean- 
ing that the threshold mass separating stars that end up 
as white dwarfs from those that end up as supcrnovae 
increases, as seen in Fig. [7] 

This may have interesting observational implications. 
In particular, in recent years there has been impor- 
tant progress identifying progenitors of observed super- 
novae. Hints are accumulating that stars of initial mass 
8M© — 9M© do, in fact, end up as supernova explo- 
sions. A particular case is supernova SN 2003gd (Type 
II-P), for which the observations fav or a red supergiant 
progenitor of mass ~ 8M© — 9 M© (jSmartt et al.ll2004l ; 
iVan Dvk et "aDf2003l . f2006rT . 

The situation with observations continues to improve 
at a rapid pace. For example, very recently a progenitor 
of supernova SN 2008bk (also a Type II-P) was identi- 
fied to be a red supe rgiant with the mass of 8.5 ±1.0 M© 
(|Mattila et al.ll2008h . Likewise, a recent study of 20 pro- 
genitors of Type II-P supcrnovae finds that the minimum 
stellar mass for these stars is 8+} 5 M© (|Smartt et al.l 
I2008T) . Referring back to Fig. we see that values 




1 1 — I — I 1 1 1 — I 1 1 1 — I 1—3* 

7 8 9 10 11 12 M/Mq 

Fig. 7. — The fate of the stars of different masses as a function 
of the neutrino magnetic moment (in units of 10 -12 ^b)- The low 
mass stars end up as carbon-oxygen white dwarfs; heavier stars end 
up as oxygen-neon-magnesium white dwarfs (some of which result 
in electron capture supernovae); finally, still more massive stars 
undergo collapse in their iron core. As the magnetic moment is 
increased, the mass thresholds between the different types change 
and a new type - a star undergoing a explosive thermonuclear 
runaway of its degenerate CO core - appears. 



\i v > 0.35 x 10~ 1c Vb are already disfavored. It seems 
highly likely that more observational progress is to fol- 
low. 

The interpretation of the observations relies, among 
many factors, on the stellar models used to infer the mass 
of the progenitor from its color and luminosity. Hence, 
further refinement of the stellar models would also be 
beneficial. To this end, it must be mentioned that the 
magnetic moment cooling by itself does not change the 
relationship between the surface characteristics and the 
mass: our simulation shows changes in the central region 
of the star at late stages of the evolution, but not in the 
surface characteristics. 

4.4. A New Type of Supernova 

The appearance of a new type of supernova as a con- 
sequence of the additional cooling is perhaps the most 
remarkable finding of our simulations. Should such ob- 
jects have distinct characteristics not observed in actual 
supernova explosions, a robust bound on the magnetic 
moment, [l v < 0.2 x 10~ 10 /ib, would result. 

We again stress that we were unable to follow the de- 
velopment of the explosion beyond the ignition of the 
degenerate CO burning. Such calculations would be 
highly desirable. Provisionally, we could consider a plau- 
sible model of a Type la supernova exploding inside a 
~ 9M© — 11 M© red supergiant. A simple estimate we 
performed shows that, especially for stars in the low 
mass end of the "CO explosions" bin, the amount of 
energy available upon burning of the entire CO core is 
enough to disrupt the star, though a detailed hydrody- 
namic study to confirm this assessment has yet to be per- 
formed. The resulting explosion would leave no neutron 
star remnant, and no neutrino signal would be observed. 
The amount of 56 Ni and 56 Co could be close to what 
is produced in a Type la explosion, ~ 0.4 M© — 0.7 M© 



8 



Hcgcr, Fricdland, Giannotti, & Cirigliano 



- 




— i — 




10 








9.9 M„. u=50 












































; — 

; 
































' 










: 
















C12 

N14 




□ -c 


















016 

Ne20 

Mg24 




















= 


Si28 




-5 






































S32 

Ar36 






I A 
















Ca40 








-\ i :\ 






— r— i"Ny 








-10 



D.D 0.5 1.0 1.5 1.7 3 6 9 














---v- 










































10.5 Mg, ^ = 


50 














































- 1 












l\ 




































































He4 




S -2 
E 


















C12 

N14 






















015 

Mg24 














1 




L 






-3 
















= 


Si28 






















Ar36 




















Ca40 






: B 




















-A 





















0.0 0.5 1.0 1.5 1.7 5 6 9 



Fig. 8.— Final structure of 9.9 Mg (Panel A) and 10.5 Mq 
(Panel B) stars for fj, v = 0.5xl0~ 10 (j,b a t t ne time of thermonu- 
clear runaway at the center. 



(|Blinnikov et al.ll2006[ ). This would significantly exceed 
the amounts ejected in a standard core collapse explo- 
sion, ~ 0.1 M - 0.15 M Q , or in po ssible ECSNe from ou r 
ONeMg WD bin. As discussed bv lWooslev et all (12002ft . 
the decay of 56 Co to 56 Fe (half-life 77 days) is obser- 
vationally very significant, giving a "radioactive tail" to 
the light curve. Thus, excessive amounts of 56 Co could 
be seen. For the current accuracy of the determina- 
tion of 56 Ni mass f rom observations see, e.g., Fig. 9 in 
ISmartt et all (j2008f) . 

Notice that in the previous argument we concentrated 
on the low-mass end of the "CO explosions" bin. In fact, 
the pre-explosion composition of the stars in the "CO 
explosions" bin varies with the mass of the star. Two 
typical examples are shown in Fig. [8] We see that only 
the inner fraction of the core has a CO composition. This 
CO inner core is surrounded by either a ONeMg layer at 
lower masses, or a silicon-sulfur (SiS) layer (from ONeMg 
shell burning) at higher masses. We find that the size of 
the central CO "nugget" decreases as the mass of the 
star increases. This will affect densities, nucleosynthesis, 
and energetics of the explosion. 

At the upper mass limit of the "CO explosions" bin, 
the energy may eventually become insufficient to dis- 
rupt the star; in this case, a burning pulse may eject 
the outer layer just before the final core collapse SN, 
making it rather bright at the peak due to the inter- 
action of the core collapse SN ejecta and the previously 
ejected layers. This mechanism is similar t o the model for 
pair- in stability supernovae presented in IWooslev et alJ 
(|2007f ). The resulting SN would have little 56 Ni to power 
the tail of the light curve. 



In a few cases at the upper mass limit we also observe 
the formation of an "iron" layer from silicon shell burn- 
ing; these cases have only a tiny (~ 0.1 M Q ) CO core left 
in the center and collapse as core collapse SNc before 
igniting the central CO nugget; they would probably be 
indistinguishable from regular CCSNe. 

The detailed information from our simulations is col- 
lected in Table [Cl] We can see that a significant fraction 
of the supernovae in the "CO explosions" bin do have 
sizable CO cores, enough to disrupt the stars and give 
the characteristic excess 56 Co observational signature. 

Lastly, we note r ecent observations of t he Kepler's su- 
pernova remnant ([Reynolds et al.l l2007f ) 3 . The abun- 
dance of iron in the shocked ejecta together with no evi- 
dence for oxygen-rich ejecta indicate that this object un- 
derwent a thermonuclear explosion. Curiously, a detailed 
analysis described in that paper suggests the possibility 
of a Type la explosion in a more massive progenitor. 
That the explosion occurred in the environment that is 
not low in metallicity is especially interesting. Clearly, 
further observational studies are very desirable. 

4.5. Other Effects 

There are several other effects of the extra cooling 
that arc worth mentioning. For example, the result- 
ing composition of the core-collapse supernova progeni- 
tor is somewhat modified. Our calculations of the 12 M Q 
star show the size of the CO core decreases by 8 % with 
\i v = 0.5 x 10~ 10 /is- This small change could be of inter- 
est, since the explosion mechanism is a sensitive function 
of the chemical composition. 

Another e ffect is the change of the presupernova neu- 
trino signal (|Odrzvwolek et aLll2004f ). This occurs be- 
cause the last stage before the iron core collapse - Si 
burning - proceeds significantly faster since the energy 
escapes the core more easily. For 12 M and /i„ = 
0.35 x 10~ 1( Vb, the Si shell burning time is only about 8 
hours, instead of 15 hours for fi u = 0. The pre-supernova 
neutrino signal at a future megaton-sized neutrino detec- 
tor could be shorter with higher luminosity. 

Other effects may be revealed upon a more detailed in- 
vestigation. For example, the extra cooling would modify 
the nucleosynthetic yields of various elements. 

5. CONCLUSION AND PERSPECTIVES 

For certain kinds of new physics, useful astrophysi- 
cal constraints can be obtained using simple order-of- 
magnitude analytical estimates. Examples are provided 
by various theories with extra dimensions (see , e.g. , 
Cullen fc Perelsteinl 119991 iFriedland fc Giannottil 120081 : 
Hannestad fc RaffeltJ 12003 ) . In contrast, the treatment 
of the problem given here requires a detailed analysis 
with a sophisticated stellar evolution code. 

As we have seen, stars with masses in the range 7 Mq < 
AI < 18 M Q are indeed sensitive to the neutrino magnetic 
moment. The additional cooling it introduces does not 
simply change the evolutionary time scales of these stars, 
but leads to the appearance of qualitatively new features 
in the evolution. The physics involved is very rich and 
more detailed investigations of many aspects are clearly 
desirable. We have undertaken here the first step to ex- 

3 We thank Carles Badenes for bringing this interesting study 
to our attention. 
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plore the subject. Let us now summarize our main results 
and consider some perspectives. 

1) We have shown that the magnetic moment of the 
size \x v ~ 5 x 10 _11 /zb, which is below the experimen- 
tal bound, would modify the evolution of stars with the 
masses < 17 M . The evolution of heavier stars is un- 
changed. This is a consequence of the fact that heavier 
stars have higher central temperatures and as a result 
"miss" the region of the parameter space where the ex- 
tra cooling operates. 

2) The minimal value of \x v to which the star is sen- 
sitive depends on the initial mass. For stars whose ini- 
tial mass is about 9 — 10 M Q , this value is as low as 
2 x 10~ n /iB- The sensitivity of massive stars to the 
neutrino magne tic moment is thus comparable to that 
of the HB stars (jRaffelt et al.lll989D . Notice that iRaffeld 
(|1990f ) found that the analysis of RG cooling yields a 
more stringent bound, [i v < 3 x 10 _12 ^b- 

3) The extra cooling due to the neutrino magnetic mo- 
ment shifts the threshold masses for the different final 
stages of the evolution. The number of carbon-oxygen 
white dwarfs increases, while the corresponding number 
of the oxygen-neon-magnesium white dwarfs and core 
collapse supcrnovac decreases. These shifts are graphi- 
cally illustrated in Fig[7l This effect could be constrained 
using observations of supernova progenitors. Such data 
is now coming in at a rapid pace. 

4) According to our code, stars with initial mass 
~ 9.3 Mq undergo a thermonuclear explosion for values 
of \l v larger than 2 x 10~ u /!b- This is a new feature, 
absent in the standard model evolution, and in principle 
observable. The initial mass region for which this hap- 
pens enlarges with larger /i„, and is depicted in Fig. 
labeled as "CO explosion". For \i v = 5 x lCP n /!B wc 
predict that roughly one in five of the stars which do 
not become WDs end their lives with this thermonuclear 
explosion, whereas the others become CCSNc. 



5) The final composition of the presupernova is a func- 
tion of \i v . This might be of interest since the explosion 
mechanism is a sensitive function of the chemical com- 
position. 

6) For stars that end their lives as CCSNe, a possible 
presupernova neutrino signal at a future megaton water 
Cherenkov neutrino detector would be modified in the 
presence of the additional cooling. The neutrino flux 
would be higher and the duration shorter. 

7) Our conclusions are based on a numerical simula- 
tion performed with the KEPLER code. A reliable pre- 
diction of the threshold masses between different kinds of 
WDs and CCSN depends on, among other ingredients, 
the treatment of convection, which varies between the 
codes. A comparison of differ ent codes (standard evo - 
lution only) was performed by iPoelarends et al.l (|2008h . 
yielding slightly different predictions of the thresholds. 
In the case of nonstandard physics, a similar compari- 
son would be very useful. Additionally, one has to keep 
in mind other uncertainties that affect stellar evolution, 
such as the initial stellar composition (e.g., iTur et all 
|2007|) and the natural distr ibution of stellar rotation 
rates (e.g.. iHeger et al.ll2000t ). 

8) Beyond this list, various other aspects of the extra 
cooling could be considered, such as the changes in the 
nucleosynthesis of various elements. 
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APPENDIX 

APPENDIX A: ENERGY LOSS DUE TO STANDARD NEUTRINO INTERACTIONS 

In this appendix we give approximate expressions for the standard neutrino cooling rates as a function of density 
and temperature, in order to give an analytical description of Fig. [3^- We restrict ourselves only to the regions relevant 
to the evolution of massive stars. 

At temperatures below <5x 10 8 K and low densities (orange regio n in the middle left in Fig. [3jt) the photo process 
dominates. Energy loss per unit mass due to this process is roughly (jPetrosian et al.lll967| ) 

e photo ~ 0.1 erg g- 1 s" 1 x T| , (Al) 

where Tg = T/10 8 K. At higher temperatures ( top turquoise region in Fig. [3^) pair production instead dominates. 
The energy loss is roug hlv dFowler fc HovldfT96l 

£pair ~ 4.6X10 11 ergg- 1 s" 1 ^ e" 118 ' 5 ^ , (A2) 

Pi 

for T < few 10 9 K and 

Spai, ~ 2.8xl0 20 crgg- 1 s- 1 ^ , (A3) 

Pi 

for T > few 10 9 K, where p 4 = p/10 4 gcm- 3 and T w = T/10 10 K. 

Both the above equations apply only for non-degenerate plasma. When degeneracy becomes important the pair 
emission is exponentially suppressed by the electron chemical potential. In this regime, the main loss mechanism is 
the neutrino emission via plasmon decay (green region in the middle of Fig. [3^)- This rate was first calculated by 
llnman fc Rudermanl |l964). For densities such that w p i < T, the energy loss due to this process is approximately 

Eplasma - 0.0127 ergg" 1 s- 1 If p\. (A4) 
At very high densities (u> p \ ;§> T), this emission is exponentially suppressed by the plasma frequency and the 
brcmsstrahlung process becomes the dominant one. 
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APPENDIX B: ENERGY LOSS DUE TO NON-STANDARD NEUTRINO INTERACTIONS 

If neutrinos have a non-vanishing magnetic moment, the production mechanisms (pair, plasma, photon, 
brcmsstrahlung) would receive an additional electromagnetic contribution from the interaction term of Eq.JT]), as 
shown in Fig. [5] We restrict ourselves only to the processes relevant for the evolution of massive stars, namely plasma 
and pair production (the photo process dominates neutrino energy loss in a region where the star cools mainly via 
radiative energy loss). 

Plasma processes 

The standard and electromagnetically induced plasmon emission rates (e p i a sma and £pi asma ) depend differently on the 
plasma frequency, and hence on the density. This behavior can be anticipated on the grounds of dimensional analysis. 
The energy loss rate can be schematically written as (decay rate) x (photon energy) x (number density of photons). 
Therefore the ratio of the standard e p i asma and the add itional nonstandard £p l asma cooling rates is equal to the 
corresponding ratio of the plasmon decay rates, T^/T (cf. iFriedland fc Giannottil [20081 ). The latter, in turn, can be 
estimated by dimensional analysis. The rates are proportional to the squares of the corresponding coupling constants, 
p v or G-p/e. The remaining dimension should be fixed by lu p \ alone. (In the rest frame of the plasmon, T cannot enter.) 
Therefore, T^/T ~ p 2 , / [G 2 F / oi]uj~ 2 . In the region of temperature and density of interest for us, a detailed calculation 
gives (lHaft et al.lll994f ) 

£ plasma nolo /10keV\ 2 / D1 \ 

~ 0.318 p 12 , (Bl) 

£ plasma \ ^pl / 

where the plasma frequency is numerically given by, e.o.. lRaffeltl (fl996h . 

** = 28 - 7cV [l + (1.019 Tifw ' (B2) 

with p in gem -3 . Therefore at low densities, p <C 10 6 gcm~ 3 , the ratio between the electromagnetic induced plasmon 
emission rate and the standard one is inversely proportional to the density, whereas at high density p 3> 10 6 gcm -3 , 
it is inversely proportional to p 2 / 3 . Finally, taking for example p„/pb = 10~ 10 , equation (|B1[) shows that, for density 
sufficiently low, the non-standard plasma process can be much larger than the standard one. Therefore, when the 
non-standard neutrino interactions are turned on, the region where the plasma process is dominant becomes larger, 
especially at low densities (see panel (b) and (c) of Fig. [3]). 

Pair production 

The energy loss per unit time and volume in neutrino pair production is given by the transition probability for 
e + e~ — > vv annihilation multiplied by the neutrino pair energy and integrated over the density of both initial and 
final states. This reads 

= (2^ / e(^ + le(W + l ^ + E ^ ^ ' (B3) 

where p is the chemical potential, u ro i is the electron-positron relative velocity and cr pair represents the spin averaged 
cross-section for e + e~ — > vv annihilation. The quantity i?i£ , 2freic pair is Lorentz invariant and is given by (j>i = (Ei, p i ), 
q t = (Wi,qj)) 

Ek^O* = \ J (^||^|M| 2 (2*)^(Px +P 2 - ft - fc) (B4) 

in terms of the spin- averaged squared invariant amplitude \M\ 2 . If the pair production proceeds via neutrino magnetic 
moment, we find (s = (p± + P2) 2 , t = (p± — qi) 2 , u = (pi — q2) 2 ) 

\M,J 2 = ^[s 2 -(u-t) 2 } , (B5) 



where p 2 represents the effective dipole moment: 



Integrating over the invariant phase space leads to 

^pair = < 

12 

or equivalently WroicP^ 11 = (ap w ) 2 /3 (1 + 2m 2 /s). Performing the angular integrations in Eq. IB4I we then find 

OP a ir _ «M 2 , r , F f]F n n (E 1 +E 2 )(2m 2 c +E 1 E 2 ) 
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Fig. B9. — Contours of the ratio Qf^ 11 / of non-standard over standard energy loss via neutrino pair production as a function of 
temperature and density. The effective neutrino magnetic moment is fixed at /zio = 1. 



showing tha t the expression for t he luminosity factorizes in the product of two one-dimensional integrals. 
Following IBeaudet et all (|1967f) . we now define the dimcnsionlcss variables 

A = T/m c , v = n/T , (B9) 

and the dimensionless integrals 

roo 2n+l( 2 _ \-2\l/2 

G±(A,i/)=A 3 + 2 " / dx- } + ,\> ■ (BIO) 

The luminosity due to non-standard pair production reads 



)P air = a ^ m c 
6tt 4 



2 ™7 . 

2G_ 1/ , 2 G + 2G G_ 1 y 2 + G G 1 j 2 + G 1 ^ 2 G , (Bll 



and is a function of the temperature and the chemical potential. In order to have a function of temperature and 
density, the chemi cal potential must be solved in terms of temperature and density. This requires the solution of the 
following equation(Bcaudct et al.lll967| ) 

p = 2.922 x 10 6 gcm- 3 (G~ - G+) , (B12) 

which cannot i n gen eral be inverted analytically in terms of p. The numerical solutions are plotted in Fig. 1 of 
IBeaudet et ail (jl967l ). 

In order to assess the im pact of non-standard cooling, the above result has to be compared to the Standard Model 
emission rate (IDicuslll972D : 

G 2 m 9 r 

Qsm ( v e) = iftZi'l^^'v ~ 2G A )[G ^-1/2 + G_ 1 / 2 G ] + 9C V [G 1 ^ 2 G Q + G G 1/2 ] 

+ (Cy + Ca)[^1 @t/2 + ^1/2^1 _ ^1 ^-1/2 _ ^1/2^0 ~ ^0 @t/2 ~ ^-l/a^ll } ■ (B13) 

Here Cy = 1/2 + 2sin 2 6'w and Ca = 1/2, with the weak mixing angle numerically given by sin 2 0w = 0.23. The 
energy loss due to Vn and v T is obtained form Qg^J (f e ) by replacing Cy — > Gy — 1 and Ga — > —Ca- 

A simple comparison of the standard and non-standard pair production amplitudes suggests that the two contri- 
butions are of comparable size if p v ~ 10~ 10 ^b- Using the full expressions of Eqs. [Bill and IB13I we plot in Fig. IB9I 
the ratio Q5 alr /QsM f° r = 10~ 10 Mb in the region of interest in the p — T plane. The results show that indeed 
the non-standard loss never exceeds 30% of the electro-weak one, being most important in the non-relativistic and 
non-degenerate region. From the behavior of G^(A, v) in the relativistic and/or degenerate limit one can simply verify 
that Q^ air /QgM r — > in these physical regimes, as shown in Fig. IB9I 

For the numerical implementation in the stellar evolution code we have used the simple fitting formula 



QjT r 

n pair 



(p, T) = p( Q x \n(p) + o(p)(]Qg 10 T - 7.5) 2 + 6(p)(log 10 T - 7.5) io j , (B14) 

n(p) = 2.959105 + 0.0021515(log 10 p - 3) 5 , (B15) 

a(p) = 0.362925 + 0.000408304 p - 436811 , (B16) 

b(p) = 0.00998357 + 2.76174 • 10~ 7 p - 5018889 , (B17) 
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where piQ is defied by p v = /xio 10 10 pb, 2 1 = T/K, and p = (p//x e )/(gcm 3 ). The fitting formula is accurate at the 
level of 5% or better in the region defined by 10 7 5 < T < 10 10 and 10 3 < p < 10 10 . 

Finally, in the region where the pair process dominates neutrino emission, the additional energy loss per unit time 
and mass is approximately given by: 

< air ~ 1.6 x 10 11 crgg" 1 s" 1 x ^ e -n8.5/T 8 t (B18) 

Pi 

where pio = p v /10 10 pb- Comparing with Eqs. (|A2j) and (|A3|) we see that this energy loss is never more that ~ 30% 
of the standard model value, if p v is below the experimental bound. 

APPENDIX C: DETAILED LISTING OF MODEL RESULTS 
Table [CTl summarizes our calculations and results for different values of p v . 
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Neutrino Magnetic Moment in Massive Stars 

Table CI 
Overview of Model Results 



M/Mq Comment 

flu = 0.00x10" 10 

AGB, CO WD 

AGB, CO WD with 0.40 Mq NeMgO outer layer, 0.70 M s CO center 
AGB, ONeMg WD, 3rd dredge-up, vast C shell burning, 0.02 Mq CO core left in center 
AGB, OMeMg WD, 3rd dredge-up 
3rd dredge-up, may make ECSN 

NeO thermonuclear runaway, weak explosion; fallback could still be possible 
CCSN, 0.21 Mq ONeMg core, explosive ignition inside Fe shell 
CCSN 



Hv = 0.10x10- 1[ Vb 

AGB, CO WD 

AGB, ONeMg WD, 3rd dredge-up, vast C shell burning, up to 0.18 Mq CO core left in center 
AGB, ONeMg WD (3rd dredge-up) 
ONeMg WD or ECSN 

CCSN, ONeMg core with SI Shell burning, 

CCSN, central NeO ignition, weak explosion, fallback of Fe core 
CCSN 



ii v = 0.20x10- 11 Vb 

AGB, CO WD 

AGB, ONeMg WD 

O/Ne/C thermonuclear runaway 

0.32 M SiS core (10% O), CCSN 

Si explosion, 3rd dredge-up of He shell, CCSN 

CCSN, Si detonation 

CCSN 

CCSN, 0.20 Mq SiS nugget at center, should make Si detonation 
CCSN 



ii v = 0.35x10- 1( Vb 

AGB, CO WD 

AGB, 1.08 - 1.13 M Q CO WD with 0.10 - 0.40M Q ONeMg layer 
AGB, 1.15 - 1.32 Mq ONeMg WD with ~ 25% C in center 
AGB, 1.35 Mq ONeMg, 25% C in center, deep 3rd dredge-up 

C thermonuclear runaway, 0.63 — 0.50 Mq CO nugget, inside 1.37 — 1.00 Mq ONeMg core 

C thermonuclear runaway, 0.25 — 0.03 Mq CO nugget in center 

CCSN 



ii v = 0.50x10- 1c Vb 





4 


.0 


AGB, CO WD, 0.2% He left close to center 




4.5 


- 7 


.5 


AGB, CO WD 




8.0 


- 8 


.3 


AGB, CO WD, starting and increasing C shell burning 




8.4 


- 9 





AGB, CO WD, 0.20 - 0.60 M Q layer with Ne > C 




9.1 


- 9 


.6 


AGB, ONeMg WD, 0.60 - 0.49 M CO core 






9 


7 


AGB, ONeMg WD, 0.52 Mq CO core, no 3rd dredge-up 






9 


.8 


C thermonuclear runaway, O.7M CO core, 1.36 Mq He core, 3rd dredge-up, NeO shell 


9.9 - 


- 10 


.4 


C thermonuclear runaway, 0.60 — 0.49 Mq CO core, 2.46 - 


- 2.62 M He core, NeO shell 


10.5 - 


- 10 


.9 


C thermonuclear runaway, 0.30 — 0.10 Mq CO core, 2.66 - 


- 2.78 M He core, SiS shell 




11 


.0 


C thermonuclear runaway, 0.02 Mq CO core, 2.810 M He core, Fe shell =>• CCSN 




11 


.1 


C thermonuclear runaway, 0.05 Mq CO core, 2.80 Mq He 


core, SiS shell 




11 


2 


C thermonuclear runaway, 0.05 Mq CO core, 2.82 Mq He 


core, Fe shell =► CCSN 




11 


.3 


C thermonuclear runaway, 0.17Mq CO core, 2.76 Mq He 


core, SiS shell 




11 


.4 


C thermonuclear runaway, 0.35 Mq CO core, 2.62 Mq He 


core, NeO shell 




11 


.5 


C thermonuclear runaway, 0.15 Mq CO core, 2.74 Mq He 


core, SiS shell 


11.6 - 


- 25 


.0 


CCSN 





4.0 - 7.1 
7.2 
7.3 

7.4 - 8.9 
9.0 - 9.1 
9.2 
9.3 

9.4 - 25.0 



4.0 - 7.4 
7.5 - 8.9 
9 - 9.1 
9.2 
9.3 
9.4 

9.5 - 25.0 



4.0 - 7.7 
7.8 - 9.2 
9.3 
9.4 
9.5 
9.6,9.8 
9.7,9.9 
10.0 
10.5 - 25.0 



4.0 - 8.0 

8.1 - 8.4 

8.5 - 9.4 

9.5 

9.6 - 9.9 
10.0 - 10.4 
10.5 - 25.0 



